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Tests of hard and soft QCD with e“^e Annihilation Data 
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Experimental tests of QCD predictions for event shape distributions combining contributions from hard and soft 
processes are discussed. The hard processes are predicted by perturbative QCD calculations. The soft processes 
cannot be calculated directly using perturbative QCD, they are treated by a power correction model based on 
the analysis of infrared renormalons. Furthermore, an analysis of the gauge structure of QCD is presented using 
fits of the colour factors within the same combined QCD predictions. 


1. INTRODUCTION 

Hadron production in e''"e“ annihilation is a 
well-suited process for the study of the interplay 
between hard and soft QCD. The lack of inter¬ 
ference with the initial state improves the under¬ 
standing of the hadronic final state in terms of 
QCD. The availability of data sets at many en¬ 
ergy points with centre of mass (cms) energies y/s 
between 14 and 209 GeV allows to separate the 
hard and soft contributions. 

The hard QCD processes, i.e. the radiation 
of gluons from the quarks produced at the elec- 
troweak vertex and subsequent higher order pro¬ 
cesses are well understood in perturbative QCD 
(pQCD) i. The production of multi-jet events 
is clear evidence for the existence of gluons and 
is successfully described by pQCD. 

The soft contributions to hadron production in 
e"'"e“ annihilation are connected with the transi¬ 
tion from the partons (quarks and gluons) gener¬ 
ated by hard processes to the hadrons observed 
in the experiment. This transition is referred 
to as hadronisation and takes place at energy 
scales of approximately a light hadron mass, i.e. 
0(100) MeV, where pQCD calculations are not 
reliable due to the rapid growth of the strong 
coupling as- The soft contributions are known 
to scale like an inverse power of the cms energy; 
for most observables the scaling is l/\/s. 


2. POWER CORRECTIONS AND IN¬ 
FRARED RENORMALONS 

At energy scales Q = Aqcd perturbative evolu¬ 
tion of as (<5) breaks down completely due to the 
Landau pole at Aqcd- It is therefore impossible 
to attempt pQCD calculations of soft processes 
without further assumptions. 

In order to clarify which assumptions are 
needed a phenomenological model of hadronisa¬ 
tion is considered. The longitudinal phase space 
model or tube model goes back to Feynman (P. 
One considers a qq system produced e.g. in e+e" 
annihilation in the cms system. The two primary 
partons move apart with velocity u/c ~ 1. In 
such a situation the production of soft gluons will 
be approximately independent of their (pseudo-) 
rapidity rj' = — log(tan(0i/2)), see figure ^ 

The change to e.g. the observable Thrust 1 — T 
due to the production of a soft gluon at angle 0i 
with transverse momentum kt^i is A(1 — T)i ~ 
■ This observation is generalised to 
write the soft or non-perturbative contribution to 
the mean value of the 1—T distribution as follows: 

(1 - T)np = J ■ J (1) 

The function <i)(fct) ^ as(fct) is the distribu¬ 
tion of soft particles in fct- The first integral in 
equation (H) is summarised as ao jQ independent 
of the observable and the second as a constant 
ci-T dependent on the observable. The quan¬ 
tity ao ^ / as(fct)dA:t can only exist when as(fct) 
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Figure 1. Sketch of a qq system in rj' — kt space. 


is identified with a non-perturbative strong cou¬ 
pling which is assumed to be finite at low 

kt around and below the Landau pole. 

A more formal approach to the origin of soft 
contributions is the study of infrared renor- 
malons, i.e. the divergence of asymptotic pQCD 
predictions due to the integration of low momenta 
in quark loops in gluon lines Q. The infrared 
renormalon divergence of the 0{ag) term in a 
pQCD prediction is factorial in n: ^ 

(2/3o/p)"n!Q;g^^. With Sterlings formula to re¬ 
place n! one finds ^ (2/3oQ!s/p)"ri”e“"’as- 

The convergence of the series is optimal for n = 
p/(2/3oQ;s). Based on this relation one finds 
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where the first order relation between as and 
Aqcd has been used. The result shows that in¬ 
frared renormalon contributions to pQCD predic¬ 
tions scale like Q~p similar to the soft contribu¬ 
tion studied in the tube model, see equation ([|). 

The power correction model of Dokshitzer, 
Marchesini and Webber (DMW) extracts the 
structure of power correction terms from anal¬ 
ysis of infrared renormalon contributions ^]. 
The model assumes that a non-perturbative 
strong coupling exists around and below the 


Landau pole and that the quantity ao(/ii) = 
l/^i fg' ag^(kt)dkt can be defined. The value of 
/ri is chosen to be safely within the perturbative 
region, usually /ri = 2 GeV. 

The main result for the effects of power cor¬ 
rections on distributions A(y) of the event shape 
observables 1 — T, Mh and C is that the pertur¬ 
bative prediction Fpt( 2/) is shifted [^|-0: 

F{y)=FMy-CyP) (3) 

where Cy is an observable dependent constant and 
P ~ Mfii/Q{ao{y,i) — as) is universal, i.e. in¬ 
dependent of the observable |^. The factor P 
contains the 1/Q scaling and the so-called Milan- 
factor M which takes two-loop effects into ac¬ 
count. The non-perturbative parameter ao is 
explicitly matched with the perturbative strong 
coupling as. For the event shape observables 
and Bw the predictions are more involved and 
the shape of the pQCD prediction is modified in 
addition to the shift Q] . For mean values of 1 — T, 
Mh and C the prediction is: 

{y) — (2/)pt P CyP (4) 

For (Bt) and (Bw) the predictions are also more 
involved. 


3. EXPERIMENTAL TESTS 


Distributions and mean values of event shape 
observables measured in e“'"e“ annihilation at 
many cms energy points between 14 and more 
than 200 GeV have been used recently by sev¬ 
eral groups to study power corrections M- 
This report will present mainly results from^. 
In all of these studies pQCD predictions in 
C>(ag)-|-NLLA |p^-p^ are used for event shape 
distributions while 0{ag) |^| predictions are em¬ 
ployed for mean values. The pQCD predictions 
are added together with the power correction pre¬ 
dictions. The resulting expressions are functions 
of the strong coupling as(Mz) and of the non- 
perturbative parameter ag. The complete pre¬ 
dictions are compared with data for event shape 
distributions or mean values as published by the 
experiments. 

Figure || shows the results fitting the data 
for 1 — T distributions measured at cms ener- 
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1/a da/d(l-T) 



1-T 


OPAL 189 GeV 
L3 189 GeV 
OPAL 183 GeV 
DELPHI 183 GeV 
L3 183 GeV 
OPAL 172 GeV 
DELPHI 172 GeV 
L3 172 GeV 
OPAL 161 GeV 
DELPHI 161 GeV 
L3 161 GeV 
OPAL 133 GeV 
ALEPH 133 GeV 
DELPHI 133 GeV 
L3 133 GeV 
OPAL 91 GeV 
ALEPH 91 GeV 
DELPHI 91 GeV 
L3 91 GeV 
SLD 91 GeV 
AMY 55 GeV 
JADE 44 GeV 
TASSO 44 GeV 
JADE 35 GeV 
TASSO 35 GeV 
MARK2 29 GeV 
HRS 29 GeV 
TASSO 22 GeV 
TASSO 14 GeV 


Figure 2. Scaled distributions of l — T. The solid 
lines show the result of the fit of the combined 
pQCD and power correction prediction. The dot¬ 
ted lines represent an extrapolation of the re¬ 
sult i. 


gies from 14 to 189 GeV. The data from exper¬ 
iments at ^/s < Mz are corrected for the effects 
of e+e" ^ bb events [^. The data in the fitted 
regions are well reproduced by the theory (solid 
lines). The dotted lines present extrapolations 
outside of the htted regions using the fit results 
for as{Mz) and ao which describe the data rea¬ 
sonably well. The agreement between theory and 
experiment is similar for the other observables. 

Figure ^ (solid line) presents the result of a fit 
to measurements of (M^) at y/s = 14 to 189 GeV. 
The data are well described by the htted theory. 



Figure 3. The cms energy dependence of (1 — T) 
is shown. The solid lines show the result of the ht 
of combined 0(a|) QGD calculations with power 
corrections while the dashed lines indicate the 
perturbative contribution |^. 


The dotted line in hgure ^ shows the perturbative 
part of the prediction; one observes that the soft 
(power correction) contributions increase at low 

yi. 

The combined results from the analysis of 
distributions are as{Mz) = O.lllllo;™^® and 
ao(2 GeV) = 0.579to;o72- The relatively small 
value of as{Mz) compared e.g. to the world av¬ 
erage as{Mz) — 0.1184 ± 0.0034 |]^ is due to 
small results from Mh and B^^. The individual 
and the combined results are shown as points with 
one-standard-deviation error ellipses in figure ^ 
a). 

From the analysis of mean values the com¬ 
bined results are as{Mz) = 0.1187taoo2o 
ao(2 GeV) = 0.4851 o q 45 in reasonable agreement 
with the results from distributions. Figure ^ b) 
presents the individual and combined results from 
hts to mean values as points with one-standard- 
deviation error ellipses. Gombining both analyses 
yields the hnal results: as{Mz) = 0.11711q;qq2o 
and ao(2 GeV) = 0.513;°;°^® . 

Figure ^ shows a summary of results for 
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Distributions 



as(Mz) 


Mean Values 





Figure 4. Results for as{Mz) and ao{2 GeV) 
for distributions (top) and mean values (bottom). 
The error ellipses show one standard deviation 
uncertainties (38% CL) 


cts{Mz) and ao(2 GeV) obtained in compara¬ 
ble analyses by several groups mi- One ob¬ 
serves reasonable agreement between the individ¬ 
ual measurements. The comparable results for 
Mh and i?w from Q and 0 differ, because 
in hadron mass effects on power corrections 
are considered |^. The non-perturbative param¬ 
eter ao appears universal within about 20% as 
expected while the values for as{Mz) are com- 





Figure 5. The figure shows a summary of results 
for as{Mz) and ao(2 GeV) from power correction 
studies (|-0. The results labelled PPTSL are 
from 1^. 


patible with the world average @- 

4. QCD COLOUR FACTORS 

The analysis of event shape distributions de¬ 
scribed above can be generalised to study the 
gauge structure of QCD |^]. The C>(a|)-|-NLLA 
pQCD predictions can be decomposed into ad¬ 
ditive terms which are proportional to products 
of the QCD colour factors. This analysis consti¬ 
tutes an experimental test of the gauge symme¬ 
try of QCD via radiative corrections. The values 
of the colour factors correspond to the relative 
size of contributions from three fundamental pro¬ 
cesses i) gluon radiation from a quark (Cp), ii) 
conversion of a gluon in a qq pair [T-pUf) and iii) 
radiation of a gluon from a gluon (triple gluon 
vertex TGV) (Ca). Higher order processes like 
the quartic gluon vertex are not accessible with 
the available pQCD calculations. 

The values of the colour factors are determined 
in the theory when a particular gauge symmetry 
group is chosen under which the Lagrangian of 
theory must be locally invariant. The local gauge 
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Figure 6. The figure at the top shows results for 
as{Mz), ao(2 GeV) and one of the colour factors 
for individual observables as indicated. The error 
bars show total uncertainties and the dotted lines 
show the expectation from SU(3) QCD. The fig¬ 
ure at the bottom presents combined results for 
Ca and Cp from fits to as{Mz), Ca and Cp with 
1 — T and C. The square and triangle symbols 
indicate expectations for Ca and Cp for several 
symmetry groups 0 - 


symmetry for QCD is SU(3) with Cf = 4/3, 
Ca = 3 and Tp = 1/2. The number of ac¬ 
tive quark flavours is rif = 5 at the cms ener¬ 


gies of the currently available data. The per¬ 
turbative evolution (running) of the strong cou¬ 
pling is also sensitive to the colour factors since 
quark loops and gluon loops have the opposite 
effect. In C(as) the evolution of the strong 
coupling between two scales Q and /r is given 
by as{Q) = q;s(m)/(1 + 2/3oas(M) ln(Q/M)) with 
/3o = (IICa - 2nf)/(12^). 

The power correction calculations are also func¬ 
tions of the QCD colour factors. In the analysis 
the colour factors in the pQCD predictions as well 
as in the power corrections are considered as addi¬ 
tional free parameters to be determined from the 
data. In this way the dependence on the hadroni- 
sation model is reduced compared to traditional 
analyses relying on Monte Carlo models of hadro- 
nisation, because in the Monte Carlo models the 
colour factors are fixed [l^-p0[. 

Figure || (top) shows the results of simultane¬ 
ous fits of as{Mz), ao(2 GeV) and one colour 
factor where the other colour factors have been 
kept fixed at their SU(3) values. The fits are sta¬ 
ble with the observables 1 — T and C. Some fits 
are not stable with the observables B'p and i?w- 
In alternative fits the non-perturbative parame¬ 
ter oo was kept fixed at a previously measured 
value and is varied by its total error as a system¬ 
atic uncertainty. The results and the uncertain¬ 
ties are consistent between the two types of fit. 
The variation of ag in the fits avoids a possible 
bias which might be present when oq is fixed to a 
value measured with colour factors fixed to their 
standard values. Gombined results for the colour 
factors measured individually arerif = 5.64±1.35, 
Ca = 2.88±0.27and Cp = I.45±0.27. The corre¬ 
sponding results for as{Mz) and oq are consistent 
with previous measurements. 

Figure (bottom) presents the combined re¬ 
sults of simultaneous fits to 1 — T and C of 
as(Mz), Ca and Cp with oq fixed to previously 
measured value. Such fits with ao as a free pa¬ 
rameter turn out to be unstable mainly due to 
limited precision of the data. The results are 
Ca = 2.84 ± 0.24 and Cp = 1.29 ± 0.18 with 
a consistent value of as{Mz) = 0.119±0.010 and 
consistent with the expectation from SU(3). Sev¬ 
eral other possible gauge groups are excluded. 
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5. SUMMARY 

The analyses of power corrections discussed in 
section ^ show that it is possible to describe event 
shape distributions and mean values in e“'"e“ an¬ 
nihilation data with combined pQCD and power 
correction predictions without the need for Monte 
Carlo based hadronisation corrections. Global 
fits of the combined theory with only as{Mz) 
and ao(2 GeV) as free parameters to event shape 
data measured at = 14 GeV up to the highest 
LEP 2 energies lead to a satisfactory agreement 
with the data. The non-perturbative parameter 
ao(2 GeV) is found to be universal, i.e. inde¬ 
pendent of the observable, within the expected 
theoretical uncertainty of about 20%. 

A generalisation of the power correction anal¬ 
ysis to also vary the QCD colour factors pro¬ 
vides a measurement of the colour factors based 
on the QCD radiative corrections in the pQCD 
and the power correction predictions. The re¬ 
sults for the colour factors are consistent with the 
expectations from SU(3) as the gauge group of 
QCD while the total uncertainties are competi¬ 
tive with analyses of angular correlations in 4-jet 
final states in e“'"e“ annihilation. 

One of the legacies of LEP, SLC and their 
predecessors TRISTAN, PETRA and PEP is 
the wealth of measurements using hadronic final 
states in e“''e“ annihilation. These data made 
is possible to test and verify combined predic¬ 
tions for hard (pQCD) and soft (power correc¬ 
tions) processes. Based on these tests our under¬ 
standing of the hard and soft processes and their 
interplay has improved significantly. 

The author would like to take the opportunity 
to thank the organisers for a stimulating and en¬ 
tertaining meeting. 
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